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ABSTRACT 

The intrinsic alignments of galaxies are recognised as a contaminant to weak gravita¬ 
tional leasing measurements. In this work, we study the alignment of galaxy shapes 
and spins at low redshift (z ^ 0.5) in Horizon-AGN, an adaptive-mesh-refinement hy¬ 
drodynamical cosmological simulation box of 100 h~^ Mpc a side with AGN feedback 
implementation. We find that spheroidal galaxies in the simulation show a tendency 
to be aligned radially towards over-densities in the dark matter density field and other 
spheroidals. This trend is in agreement with observations, but the amplitude of the 
signal depends strongly on how shapes are measured and how galaxies are selected 
in the simulation. Disc galaxies show a tendency to be oriented tangentially around 
spheroidals in three-dimensions. While this signal seems suppressed in projection, this 
does not guarantee that disc alignments can be safely ignored in future weak lensing 
surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement 
with observations and other simulation works, but we find less alignment for lower 
luminosity populations. We also characterize the systematics of galaxy shapes in the 
simulation and show that they can be safely neglected when measuring the correlation 
of the density field and galaxy ellipticities. 

Key words: cosmology: theory — gravitational lensing: weak - large-scale structure 
of Universe — methods: numerical 


1 INTRODUCTION 

There is mounting observational evidence that galaxies 
are subject to ‘intrinsic alignments’, i.e. correlations 
of their shapes across large separations due to tidal 
effects that act to align them in preferential direc¬ 
tions with respect to one another (Brown et al. 2002; 
Aubert, Pichon & Colombi 2004; Mandelbaum et al. 2006a; 
Hirata et al. 2007; Joachimi et al. 2011; Heymans et al. 
2013; Singh, Mandelbaum & More 2015). These alignments 
have been identified as an important systematic to weak 
lensing measurements, with the potential to undermine 
its capabilities as a probe of precision cosmology if unac¬ 
counted for (Hirata & Seljak 2004, 2010; Bridle & King 
2007; Kirk et al. 2012; Krause, Eifler & Blazek 2015). 
To fully extract cosmological information from weak 
lensing, the alignment signal needs to be mitigated 
or marginalized over (Zhang 2010; Joachimi & Bridle 
2010; Joachimi & Schneider 2010; Troxel & Ishak 
2012). On the other hand, there is cosmological in¬ 
formation to be extracted from the alignment signal 


(Chisari & Dvorkin 2013; Chisari, Dvorkin & Schmidt 
2014; Schmidt, Chisari & Dvorkin 2015). For these reasons, 
this field is emerging as an interesting avenue for improving 
our understanding of galaxy formation and the evolution of 
the large-scale structure of the Universe. For a set of com¬ 
prehensive reviews on the topic of intrinsic alignments, see 
Troxel & Ishak (2015), Joachimi et al. (2015), Kirk et al. 
(2015) and Kiessling et al. (2015). 

The large-scale alignment signal of the shapes of lu¬ 
minous red galaxies (LRGs) with the matter density field 
has been considered the dominant contribution to low red- 
shift intrinsic alignments and the main source of con¬ 
cern for weak lensing surveys. This signal has also been 
measured in numerical simulations of galaxy formation 
(Tenneti et al. 2014, 2015). A theoretical model developed 
by Catelan, Kamionkowski & Blandford (2001) suggested 
that elliptical galaxies, supported by random stellar mo¬ 
tions, can be subject to large-scale tides, which stretch them 
along the direction of the tidal field. Alignments of LRGs 
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are well described by this model (Blazek, McQuinn & Seljak 

2011 ). 

Disc galaxies, which have significant angu¬ 
lar momentum, can suffer a torque from the sur¬ 
rounding tidal field, which correlates their orien¬ 
tations (Catelan, Kamionkowski & Blandford 2001; 
Hui & Zhang 2002; Schaefer 2009; Schaefer & Merkel 
2015). Hydrodynamical and Af-body cosmological sim¬ 
ulations have reported alignments of the spins of 
dark matter halos and galaxies with the cosmic web 
(Bailin & Steinmetz 2005; Aragon-Calvo et al. 2007; 
Hahn et al. 2007; Sousbie et al. 2008; Zhang et al. 
2009; Hahn, Teyssier & Carollo 2010; Codis et al. 2012; 
Libeskind et al. 2012; Dubois et al. 2014) and there 
is indeed observational evidence of galaxy spin align¬ 
ments (Pen, Lee & Seljak 2000; Paz, Stasyszyn & Padilla 
2008; Jones, van de Weygaert & Aragon-Calvo 2010; 
Andrae & Jahnke 2011; Tempel & Libeskind 2013). Indi¬ 
rect evidence also comes from quasar polarizations, which 
have also been found to be aligned with the surrounding 
large-scale structure (Hutsemekers et al. 2014). However, 
at the same time, some works have claimed no alignment of 
discs (Slosar & White 2009: Hung & Ebeling 2012). 

Disc galaxies have higher star formation than el¬ 
lipticals, hence color information is typically used as a 
proxy to distinguish between types and alignment mech¬ 
anisms. Despite evidence of their spin alignments, the 
shape alignment of blue (disc) galaxies has been ob¬ 
served to be consistent with null (Mandelbaum et al. 
2011; Heymans et al. 2013). Recently, results from an 
analysis of the Horizon-AGN hydrodynamical simulation^ 
(Codis et al. 2015) have suggested that this correlation of 
the shapes of blue galaxies could be present at high signif¬ 
icance at redshift of 2 = 1.2, and it could thus become an 
additional contaminant to gravitational lensing. The main 
hypothesis behind that investigation was the use of spins as 
a proxy for galaxy shapes. 

To elucidate the discrepancy between the non-detection 
in low-redshift observations of blue galaxy shape alignments 
and the prediction of significant spin alignments in simu¬ 
lations at high redshift, it is necessary to explore the re¬ 
lation between projected shapes and spins in the simula¬ 
tion in the redshift range where observational constraints 
are available. We therefore propose in this work to concen¬ 
trate on the relationship between shapes and spins in the 
Horizon-AGN simulation (Dubois et al. 2014). In particu¬ 
lar, we first explore the correlations between galaxy positions 
and orientations to gain insights into the physical processes 
that give rise to intrinsic alignments. We further measure 
projected correlations of galaxy positions and galaxy shapes, 
and auto-correlations of galaxy shapes, in a way that mim¬ 
ics the construction of intrinsic alignment correlations from 
observations. Note that all results are obtained at 2 = 0.5; 
the redshift evolution of this effect is left for future work. 

In Section 2, we describe the hydrodynamical cosmo¬ 
logical simulation used in this work to study the shape and 
spin alignments of galaxies, including the identification of 
galaxies in the simulation and the methods used to obtain 
their rest-frame colours. Section 3 discusses how shapes and 
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spins are measured from the simulation and their conver¬ 
gence properties. In Section 4, we present the estimators 
used for correlation functions between spins, positions and 
shapes used in this work. We discuss both three-dimensional 
and two-dimensional (projected) correlations. While three- 
dimensional correlations give us an insight into the nature 
of the alignment process, the ultimate observables in photo¬ 
metric weak lensing surveys are the projected galaxy counts- 
ellipticity correlations. Section 5 brings together the results 
from the simulation and the theoretical modelling of the 
alignment signal. We discuss these results in relation to 
works in Section 6 and we conclude in Section 7. 

Compared to previous work by Codis et al. (2015), in 
this work we lift the assumption that galaxy spins are a 
good proxy for galaxy shapes. Instead, we present results 
on alignments both for spins and shapes classifying galax¬ 
ies by their kinematic properties. Moreover, we concentrate 
in correlations between positions and orientations, rather 
than correlations among galaxy orientations. This allows us 
to neglect potential systematics arising from spurious align¬ 
ment correlations due to numerical issues (see Appendix A). 
Tenneti et al. (2015) measured projected shape correlations 
of galaxies in the smoothed-particle-hydrodynamics simula¬ 
tion MassiveBlack-II. This work presents shape and spin cor¬ 
relations measured in the Horizon-AGN simulation, which 
has a comparable volume to MassiveBlack-H but uses a dif¬ 
ferent technique: adaptive-mesh-refinement. As the proper¬ 
ties of galaxies might be sensitive to the numerical tech¬ 
nique adopted, the comparison of alignments measured in 
Horizon-AGN and MassiveBlack-H is crucial to inform pre¬ 
dictions of intrinsic alignment contamination to future weak 
gravitational lensing surveys. 


2 THE SYNTHETIC UNIVERSE 

In this section, we describe the Horizon-AGN simulation 
(see Dubois et al. 2014 for more details) and we explain how 
galaxy properties relevant to this study are extracted. 

2.1 The Horizon-AGN simulation 

The Horizon-AGN simulation is a cosmological hydrody¬ 
namical simulation in a box of L = 100h“^Mpc a side. 
It is run adopting a standard ACDM cosmology compat¬ 
ible with the WMAP-7 cosmology (Komatsu et al. 2011), 
with total matter density = 0.272, dark energy den¬ 
sity Ha = 0.728, amplitude of the matter power spectrum 
Us = 0.81, baryon density Hb = 0.045, Hubble constant 
Ho = 70.4kms“^ Mpc“^, and Us = 0.967. There are 1024® 
dark matter (DM) particles in the box, with a resulting DM 
mass resolution of Moivi.res = 8 x 10^ Mq. 

The adaptive-mesh-refinement (AMR) code RAM¬ 
SES (Teyssier 2002) has been used to run the simulation 
with 7 levels of refinement up to A® = 1 kpc. The refine¬ 
ment scheme follows a quasi-Lagrangian criterion. A new 
refinement level is triggered if the number of DM particles 
in a cell is more than 8 or if the total baryonic mass in a cell 
is 8 times the initial DM mass resolution. 

Star formation is modelled following a Schmidt law: 
p* = e*p/tff , where p* is the star formation rate density, p is 
the gas density, e* = 0.02 (Kennicutt 1998; Krumholz & Tan 
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2007) the constant star formation efficiency, and tg the local 
free-fall time of the gas. We allow star formation wherever 
the Hydrogen gas number density exceeds no = 0.1Hcm“^ 
according to a Poisson random process (Rasera & Teyssier 
2006; Dubois & Teyssier 2008) with a stellar mass resolution 
of M*,res = poAx^ ~ 2 X 10® Mo. 

Gas cooling occurs by means of H and He cool¬ 
ing down to 10“* K with a contribution from metals 
Sutherland & Dopita (1993). Following Haardt & Madau 
(1996), heating from a uniform UV background is imple¬ 
mented after the reionization redshift Zreion = 10. Metallic- 
ity is modeled as a passive variable of the gas, which varies 
according to the injection of gas ejecta from supernovae ex¬ 
plosions and stellar winds. We model stellar feedback using a 
Salpeter (1955) initial mass function with a low-mass (high- 
mass) cut-off of 0.1 Mq (100 Mq). In particular, the mechan¬ 
ical energy from supernovae type H and stellar winds fol¬ 
lows the prescription of STARBURSt99 (Leitherer et al. 1999, 
2010), and the frequency of type la supernovae explosions is 
taken from Greggio & Renzini (1983). 

Active Galactic Nuclei (AGN) feedback is modelled 
according to Dubois et al. (2012). We use a Bondi-Hoyle- 
Lyttleton accretion rate onto black holes, given by Mbh = 
47rQ:G^M|jjP/(Cs + where Mbh is the black hole 

mass, p is the average gas density, Ca is the average sound 
speed, u is the average gas velocity relative to the black 
hole velocity, and a? is a dimensionless boost factor. This 
is given by q = (p/po)^ when p > po and a = 1 oth¬ 
erwise (Booth & Schaye 2009) in order to account for our 
inability to capture the colder and higher density regions 
of the interstellar medium. The effective accretion rate onto 
black holes is not allowed to exceed the Eddington accre¬ 
tion rate: Msdd = 47rGMBHm.p/(ercrTc), where ax is the 
Thompson cross-section, c is the speed of light, rup is the 
proton mass, and Cr is the radiative efficiency, assumed to 
be equal to Cr = 0.1 for the Shakura & Sunyaev (1973) ac¬ 
cretion onto a Schwarzschild black hole. Two different modes 
of AGN feedback are implemented: the radio mode, operat¬ 
ing when X = MBH/Mndd < 0.01, and the quasar mode, 
active otherwise (see Dubois et al. 2014 for details). 

2.2 Galaxy catalogue 

Galaxies are identified in each redshift snap¬ 
shot of Horizon-AGN using the AdaptaHOP 
finder (Aubert, Pichon & Golombi 2004), which relies 
directly on the distribution of stellar particles. Twenty 
neighbours are used to compute the local density of each 
particle. To select relevant over-densities, we adopt a local 
threshold of pt = 178 times the average total matter 
density. Note that the galaxy population does not depend 
sensitively on the exact value chosen for this threshold. 
Our specific choice reflects the fact that the average density 
of galaxies located at the centre of galaxy clusters is 
comparable to that of the dark matter. The force softening 
is of approximately 10 kpc. Below this minimum size, 
substructures are treated as irrelevant. 

Only galactic structures identified with more than 50 
star particles are included in the mock catalogues. This en¬ 
ables a clear identification of galaxies, including those in the 
process of merging. At z = 0.5, there are ~ 146 000 objects 
in the galaxy catalogue, with masses between 1.7 x 10® and 


2.3 X 10^^ Mq. The galaxy stellar masses quoted in this pa¬ 
per should be understood as the sum over all star particles 
that belong to a galaxy structure identified by AdaptaHOP. 

We compute the absolute AB magnitudes and rest- 
frame colours of the mock galaxies using single stellar pop¬ 
ulation models from Bruzual & Chariot (2003) adopting a 
Salpeter initial mass function. Each stellar particle con¬ 
tributes a flux per frequency that depends on its mass, 
metallicity and age. The sum of the contribution of all star 
particles is passed through u, g, r, and i filter bands from the 
Sloan Digital Sky Survey (SDSS, Gunn et al. 2006). Fluxes 
are expressed as rest-frame quantities (i.e. that do not take 
into account the redshifting of spectra) and, for simplicity, 
dust extinction is neglected. Once all the star particles have 
been assigned a flux in each of the colour channels, we build 
the 2D projected maps for individual galaxies (satellites are 
excised with the galaxy finder). The sum of the contribution 
of their stars yields the galaxy luminosity in a given filter 
band. 


3 SPINS AND SHAPES OF GALAXIES 

To assign a spin to the galaxies, we compute the total angu¬ 
lar momentum of the star particles which make up a given 
galactic structure relative to the centre of mass. We can 
therefore write the intrinsic angular momentum vector L or 
spin of a galaxy as 

N 

L = X , (1) 

n = l 

where n denotes each stellar particle of mass position 

and velocity relative to the center of mass of that 
galaxy. The total stellar mass of a galaxy is given by M* = 
The specific angular momentum is j = LilMt. 
For each galaxy, we also obtain its V/a, stellar rotation 
versus dispersion, that we measure from their 3D distribu¬ 
tion of velocities. V/a is a proxy for galaxy morphology: 
low values of this physical quantity indicate that a galaxy 
is more elliptical, pressure-supported by random stellar mo¬ 
tions; high values oi V/a suggest a disc-like galaxy, with 
stellar rotation predominantly on a plane. We first compute 
the total angular momentum (spin) of stars in order to de¬ 
fine a set of cylindrical spatial coordinates (r, 0, z), with the 
2 -axis oriented along the galaxy spin. The velocity of each 
individual star particle is decomposed into cylindrical com¬ 
ponents Vr, vg, Vz, and the rotational velocity of a galaxy is 
V = vg, the mean of vg of individual stars. The velocity dis¬ 
persion of each velocity component ar, ag, Oz is computed 
and used for the average velocity dispersion of the galaxy 
= {al + al + aD/J,. 

Figure 1 shows the distribution of V/a for galaxies in 
the Horizon-AGN simulation. We see an increase of the ra¬ 
tio V/a from low-mass M* ~ 10® Mq to intermediate-mass 
M* ~ 2x 10^® Mq galaxies, for which this median of the ratio 
peaks at 1. For massive galaxies, above 2x 10^® Mq, the V/a 
ratio decreases with stellar mass. Dwarf and massive galax¬ 
ies have low V/a and hence are pressure-supported, while 
intermediate-mass galaxies have wide-spread V/a values and 
are thus a hybrid population of rotation- and pressure- 
supported galaxies. 
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We also measure the inertia teirsor of a galaxy to char¬ 
acterize its three dimensional shape. This tensor is given by 




( 2 ) 


and it is then diagonalised to obtain the eigenvalues Ai 
A 2 ^ As and the corresponding unit eigenvectors ui, U 2 
and U 3 (respectively minor, intermediate and major axis of 
the ellipsoid). Analogously, we will explore the differeirces 
in measuring galaxy shapes using Equation (2) compared to 
the reduced inertia tensor, which is defined by 


M, E' 




( 3 ) 


where is the three dimensional distance for the stellar 
particle n to the center of mass of the galaxy. The reduced 
inertia tensor is a closer representation of the shape of a 
galaxy as measured for weak gravitational lensing measure¬ 
ments, where the inner (and more luminous) region is up¬ 
weighted with respect to the outskirts. For the reduced in¬ 
ertia tensor, Tenneti et al. (2014) found that the iterative 
procedure reduces the impact of the spherically symmetric 
r~ weights, yielding shapes that are not as round as for 
the non-iterative procedure. We will consider here the sim¬ 
ple and the reduced inertia tensor cases, and we expect that 
the results from applying iterative procedures to define the 
galaxy shapes (Schneider, Frenk & Cole 2012; Tenneti et al. 
2014) will lie between the two cases considered. 

Projected shapes are obtained by summing over i,j = 
1,2 in Equations (2) and (3), and the semiminor and semi¬ 
major axes are correspondingly defined by the eigenvec¬ 
tors of the projected inertia tensors. The axis ratio of the 
galaxy, q = b/a, is defined from the eigenvalues of the pro¬ 
jected inertia tensor as the ratio of the minor to major axes 
{b = V^, a — '/Xl, where Aa is the largest eigenvalue and 
Xb, the smallest). The components of the complex ellipticity, 
typically used in weak lensing measurements, are given by 

(e+,ex) = Y^p-^[cos(2<))),sin(2<)))], (4) 

where (j) is the orientation angle, -|- indicates the radial com¬ 
ponent of the ellipticity and x is the 45 deg-rotated compo¬ 
nent. Intrinsic alignments typically manifest themselves as 
a net average e+ ellipticity around over-densities. The total 
ellipticity of a galaxy is thus e = e\ -\- e\. In this work, ra¬ 
dial alignments have negative e+ and tangential alignments 
(as expected from weak gravitational lensing) correspond to 
positive e+. No net correlation of the x component with the 
matter density field is expected, and this statistic is thus 
used to test for systematics in the ellipticity measurement 
procedures. We discuss other systematics tests of ellipticity 
correlations in detail in Appendix A. 


3.1 Convergence tests 

Poorly resolved galaxies in the simulation are subject to un¬ 
certainties in their shapes and orientations. We identify two 
sources of uncertainties. First, there is the inherent variance 
in the shapes that arises from the choice of particles used 
in the shape computation. We call this measurement noise. 



log M, (solar mass) 

Figure 1. y/cr as a function of mass for the galaxies used in this 
work. The gray circles represent the median and the variance in 
10 logarithmic bins of stellar mass. Dwarf and massive galaxies 
are pressure-supported, while intermediate-mass have a mixture 
of rotation and pressure support. The vertical red dashed lines 
represent the cuts corresponding to the different mass bins as in 
Codis et al. (2015). The horizontal red dashed lines represent our 
cuts in V/a, chosen such that there is approximately the same 
number of galaxies in each V/a bin. 


fmeas, in analogy to the shape measurement from galaxy im¬ 
ages for weak lensing. Second, there is a bias in the ellipticity 
and orientation measurement associated to the resolution of 
a galaxy, (Jres. To define the minimum number of particles 
needed to obtain the shape of a galaxy, we compare these 
two uncertainties to the shape noise that arises from the 
dispersion in the intrinsic distribution of the shapes, Ce. 

The resolution bias is determined by randomly subsam¬ 
pling stellar particles in each galaxy with > 1000 particles. 
We compare the ellipticity measured from random subsam¬ 
ples of 50,100, 300 and 1000 stellar particles for those galax¬ 
ies. The results are shown in Figure 2 and they suggest that 
a minimum number of 300 particles in each galaxy has to 
be required in order to guarantee that the bias in the el¬ 
lipticity is an order of magnitude below the shape noise: 
f^res ^ 0.1 CTe- The distributions of galaxy ellipticities are 
shown in Figure 3. 

We determine cTmeas by bootstrap resampling (100 
times) the stellar particles used for defining the inertia ten¬ 
sor from the overall population of stellar particles that make 
up each galaxy. We compute this uncertainty for galaxies 
with > 50, > 100, > 300 and > 1000 particles. Figure 4 
shows the distribution of uncertainty in the shape measure¬ 
ment (compared to the rms ellipticity of the galaxy sample 
with 1000 particles) for each minimum number of particles 
considered. The solid lines correspond to the shape mea¬ 
sured using the simple inertia tensor and the dashed lines, 
to those measured using the reduced inertia tensor. The un¬ 
certainties are typically below the 10 ~® level and they are 
smaller in the case of the reduced inertia tensor. Resampling 
has a larger impact on the simple inertia tensor, since there 
are few particles in the outskirts of each galaxy, and these 
contribute equally to the shape measurement as those par¬ 
ticles in the central part. Nevertheless, (Tmeas is very small 
compared to the shape noise and to the resolution bias, and 
hence can be neglected. Notice that in our modelling we do 
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Figure 2. Resolution bias, = (leso, 100,300 — 

eiOOOp)/o'^(eiOOo)i a® a function of the ellipticity obtained 
from the 1000 particle subsample for galaxies with more than 
1000 particles. The colours indicate the number of stellar 
particles in the subsample: > 50 (black), > 100 (red) and > 300 
(blue). The solid lines represent the results for the simple inertia 
tensor and the dashed lines, for the reduced inertia tensor. 
Subsampling has a larger impact on the simple inertia tensor due 
to the fact that particles in the outskirts contribute with equal 
weights as particles in the center. 


not include surface brightness cuts, attenuation by dust, the 
effect of the atmosphere or convolution by telescope optics. 
Also, the resampling of stellar particles is carried out over all 
stellar particles, which constitute a correlated data set. All 
of these could have a significant impact in the estimation of 
(Tmeas but their modelling is outside the scope of this work. 

We also estimated the impact of subsampling and re¬ 
sampling on the measurement of the specific angular mo¬ 
mentum of the galaxies. We found that the measurement 
noise is insignificant and that the specific angular momen¬ 
tum is less sensitive to the number of particles in the sub¬ 
sample than the projected shapes. While subsampling intro¬ 
duces significant scatter in the measurement of the specific 
angular momentum, there is a strong correlation between 
this quantity measured with 50, 100, 300 and 1000 stellar 
particles. This correlation is shown for 50 versus 1000 par¬ 
ticles in Figure 5. 

We conclude that the minimum number of particles re¬ 
quired to perform spin and shape measurements is deter¬ 
mined by the requirements placed on the projected shapes. 
A minimum of 300 particles is needed for the uncertainty 
in projected shape to be at least one order of magnitude 
smaller than the shape noise. From this section onwards, we 
will work only with galaxies that have more than 300 stellar 
particles in the simulation. Notice that Tenneti et al. (2015) 
adopted a cut on 1000 stellar particles, while Velliscig et al. 
(2015) adopt similar cuts as in our work on the number of 
stellar particles. Figure 6 shows a comparison between the 
axis ratios of galaxies in the simulation using the simple (y- 
axis) and the reduced (z-axis) inertia tensors for galaxies 
with more than 300 stellar particles. The use of the reduced 
inertia tensor clearly results in rounder shapes as a result of 
up-weighting the inner regions of galaxies. 



Figure 3. Distribution of ellipticities for galaxies with > 50 
(black), > 100 (red), > 300 (blue) and > 1000 particles (green) 
for the simple (solid) and the reduced (dashed) inertia tensors. 
The shapes obtained with the reduced inertia tensor are rounder 
than those obtained with the simple inertia tensor. 



log'ioKeas/(e^(lOOO))] 

Figure 4. Distribution of measurement uncertainties in the ellip¬ 
ticities for galaxies with > 50 (black), > 100 (red), > 300 (blue) 
and > 1000 particles (green) for the simple (solid) and the reduced 
(dashed) inertia tensors. The shapes obtained with the reduced 
inertia tensor have significantly smaller uncertainties. These un¬ 
certainties are insignificant compared to the resolution bias and 
to the shape noise in the sample. 


4 CORRELATION FUNCTIONS 

4.1 Correlations in three dimensions 

Three dimensional correlations of galaxy shapes and spins 
can give interesting insights into the formation processes 
leading to alignments. In this section, we define the corre¬ 
lation functions of galaxy orientations in three dimensional 
space. We define the orientation-separation correlation as a 
function of comoving separation r, 

7?e(r) = (|f • e(x-fr)l^) - 1/3, (5) 
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Figure 5. Correlation between the logarithm of the specific an¬ 
gular momentum, j = |L|/M*, when measured with 50 (j/-axis) 
and 1000 (a;-axis) particles, for galaxies in the simulation with 
> 1000 particles. 



bja (reduced) 


Figure 6. Axis ratios for galaxies with > 50 particles obtained 
with the simple inertia tensor (y-axis) an the reduced inertia ten¬ 
sor (a;-axis). Colours indicate the number of galaxies in each bin. 
While there is a clear correlation between the estimates obtained 
with the two different methods, the reduced inertia tensor tends 
to produced rounder shapes. 


where e is the unit eigenvector of the inertia tensor pointing 
in the direction of the minor axis and f is the unit separa¬ 
tion vector. A positive correlation indicates a tendency for 
the separation vector and the minor axis of a galaxy to be 
parallel, hence for the galaxy to be elongated tangentially 
with respect to another galaxy. Notice that gravitational 
leasing produces a similar correlation for galaxy pairs sepa¬ 
rated by large distances along the line of sight. A negative 
correlation corresponds to a preferential perpendicular ori¬ 
entation of the minor axis with respect to the separation 
vector, resulting in a net radial orientation of galaxy shapes 
around other galaxies. Our treatment of galaxy orientations 
corresponds to representing the galaxy with the normal of its 
orientation plane and does not take into account the intrinsic 


dispersion in ellipticities. In other words, both spheroidals 
and discs are infinitely thin in this approach, and the plane 
that represents them points in the direction perpendicular 
to e. Alternatively, we will also measure the analogous cor¬ 
relation of separation vectors with spins, s, defined from the 
angular momentum of a galaxy, rjs{r). 

Grid-locking (studied in detail in Appendix A) mani¬ 
fests itself as a correlation between spin and shapes with the 
directions of the simulation box. As a result, spin and ellip- 
ticity auto-correlations (presented in Appendix B) can be 
contaminated by this effect, as two nearby galaxies formed 
from the same stream of infalling gas will be correlated 
among themselves and with the grid in which forces are 
evaluated and the effect can be enhanced at small distances 
because of the correlation length of the force field. On the 
contrary, we demonstrate that position-shape, position-spin 
and position-orientation correlations presented in the main 
body of this work are not affected. 


4.2 Projected correlations 

Three dimensional information of galaxy shapes and orien¬ 
tations is costly and near-term imaging surveys will be lim¬ 
ited to correlations between projected shapes in broad to¬ 
mographic bins. In this section, we describe projected statis¬ 
tics of alignments. First, we determine the orientation and 
ellipticity (Equation 4) of a galaxy from the projected in¬ 
ertia tensor. We obtain the real-space correlation function 
of galaxy positions and shapes projecting along one of the 
coordinates of the simulation box in a way that mimics the 
construction of Landy-Szalay estimators (Tandy & Szalay 
1993) in imaging surveys. 

We label the dataset of tracers of the density field by 
D, the set of galaxies with ellipticities by S+ (for the tan¬ 
gential ellipticity component, and analogously for the x 
component) and the set of random points by R. We de¬ 
fine the correlation function of galaxy positions and ellip¬ 
ticities, ^g,+(rp,n) as a function of projected separation in 
the sky, Vp, and along the line of sight, 11. This function is 
estimated from the sample of galaxies by using a modified 
Landy-Szalay estimator (Mandelbaum et al. 2011) 

(iMrr.n) - . ( 6 ) 

E ‘a- m 

(’■p.n) 

where TZ is the responsivity factor (Bernstein & Jarvis 
2002), TZ = 1 — (e^), e+j is the tangential/radial compo¬ 
nent of the ellipticity vector of galaxy j and the sum is over 
galaxy pairs in given bins of projected radius and line of 
sight distance. S+R is the sum of galaxy ellipticities around 
random points. In the previous expression, we are implicitly 
assuming that there is a single random sample, a procedure 
which is valid when the galaxy sample and the shape sam¬ 
ple coincide. For the random sample, we use a set of points 
uniformly distributed in the simulation box with 10 times 
the density of the galaxy sample. For cross-correlations of 
different samples, generalization of Equation (6) is straight¬ 
forward (e.g. Singh, Mandelbaum & More 2015). 

This correlation function is then projected along the 
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line of sight by integration between —Ifmax < 11 < 

/ Ilmax 

dnCg+(rp,n), (8) 

■nmax 

where we take Ilmax to be half the length of the simulation 
box. We are also interested in the correlation between the 
density field and the galaxy shapes, ws+{rp), for which the 
galaxy positions are replaced with the positions of randomly 
sampled DM particles in the simulation box. 

Two sources of uncertainties can be identified when 
measuring projected correlation functions of positions and 
shapes. The first one is shot noise coming from the tracers 
of the density and shear field. The second source of uncer¬ 
tainty is cosmic variance from the limited volume of the 
simulation box. Given that the box is 100 h~^ Mpc on each 
side, we limit the measurement of correlations to up to a 
quarter the length of the box size, L. However, this does 
not guarantee that we will be free from cosmic variance for 
the measured separations. The uncertainties in the projected 
correlation functions, including both shape noise and cosmic 
variance from modes within the box, can be obtained from 
jackknife resampling over the simulation box (Hirata et al. 
2004; Mandelbaum et al. 2006b). We divide the simulation 
box in cubes of L/3 a side. In each jackknife iteration we 
remove the tracers corresponding to one of the cubes and 
compute a new estimate of the correlation. This allows us 
to estimate the variance in the projected correlations in¬ 
cluding the effect of cosmic variance. By contrast, the shape 
noise-only variance underestimates the jackknife variance by 
a factor of 1 — 4 on small scales (< lh~^ Mpc) and up to 
approximately an order of magnitude on the largest scales 
probed in this work. In comparison, the jackknife accounts 
for covariance between H bins and for the finite volume of 
the simulation box. Nevertheless, due to the limited size of 
the box, it is expected that the error bars could still be un¬ 
derestimated on large scales. We have also confirmed that 
using more jackknife regions with L/4 length a side does not 
alter the main conclusions of this manuscript. 

5 RESULTS 

5.1 Relative orientations in three dimensions 

In this section, we characterize three dimensional correla¬ 
tions of galaxy positions and orientations at 2 ~ 0.5 to shed 
light into the physical processes that can lead to alignments. 
Since three dimensional shapes are costly (see Huff et al. 
2013 for a potential method), we study on-the-sky galaxy- 
ellipticity correlations in Section 5.2. 

We use all galaxies in the simulation with > 300 stel¬ 
lar particles divided into three subpopulations by their V/a, 
which is obtained as described in Section 3. This allows us 

^ Notice that for simplicity we do not include the effect of pecu¬ 
liar velocities in transforming galaxy redshifts into distance along 
the line of sight, nor do we incorporate photometric redshift un¬ 
certainties. Both effects decrease the significance of intrinsic align¬ 
ment correlations; however, in this work, we are interested in ex¬ 
tracting as much information as possible from the simulation and 
we opt to avoid diluting the alignment signal by including these 
effects in the modelling. 


to determine whether the alignment trends depend on the 
dynamical properties of the galaxies. Disc galaxies are rep¬ 
resented by large values of V/a, while spheroidals have low 
V/a. As mentioned in Section 1, different alignment mech¬ 
anisms are expected to act on these different populations 
(Catelan, Kamionkowski & Blandford 2001). Disc galaxies 
are expected to become aligned with the large-scale tidal 
field due to torques in their angular momentum, while ellip¬ 
ticals might suffer stretching or accretion along preferential 
directions determined by large-scale tides. 

The three V/a bins are constructed to host 1/3 of 
the galaxy population in each bin {V/a < 0.55, 0.55 < 
V/a < 0.79 and V/a > 0.79). The left panel of Figure 7 
shows the three-dimensional correlation function of minor 
axis, from the simple inertia tensor, and separation vector. 
Error bars correspond to the standard error on the mean 
in each bin. We are only interested in a qualitative com¬ 
parison in this section and we thus neglect cosmic vari¬ 
ance except when explicitly mentioned. We observe that 
the minor axis orientation-position correlation is only sig¬ 
nificant in the low V/a bin, suggesting that only galax¬ 
ies with spheroidal dynamics are subject to shape align¬ 
ments. The results imply that these galaxies are elongated 
pointing towards each other, in agreement with the quali¬ 
tative behaviour expected from the tidal alignment model 
(Catelan, Kamionkowski & Blandford 2001). (Notice that, 
by construction, the fraction of spheroidal galaxies at this 
redshift is 1/3.) Galaxies with V/a > 0.55 do not have corre¬ 
lated positions and shapes. Moreover, low V/a have a lower 
degree of correlation between the direction of their minor 
axis and the direction of their spin, and this is a mono¬ 
tonic function of V/a. Finally, we have split the low V/a 
population of galaxies by the mean stellar mass; we find 
that both low (log M* < 9.5, ~ 20,000 galaxies) and high 
mass (log M* > 9.5, ~ 8,000 galaxies) galaxies are sub¬ 
ject to alignments, and that high mass galaxies are more 
strongly aligned. Spin alignments are not very significant 
among these populations, although there seems to be a small 
trend for galaxies being oriented perpendicularly to the sep¬ 
aration, i.e., tangentially around other galaxies. This also 
seen in the reduced inertia tensor of low and intermediate 
V/a galaxies. 

We also divide the galaxy population into three bins of 
mass: low mass (10® < M* < 10®'® Mq, ~ 30,000 galax¬ 
ies), intermediate mass (10® ® < M* < 10^®'® Mq, ~ 40,000 
galaxies) and high mass (M* > 10^®'® M©, ~ 7, 000 galax¬ 
ies). The low mass limit is determined by the required 
threshold in the number of stellar particles. We first note 
that with this selection, both low and high mass galaxies 
tend to have low V/a, while intermediate mass galaxies have 
V/a values more consistent with those of a disc-like popula¬ 
tion, as shown in Figure 1. In Figure 8, we show the correla¬ 
tion between spin or minor axis with the separation vector 
for the three mass bins. The left and middle panels show 
the minor axis-separation correlation constructed using the 
simple (left) and reduced (middle) inertia tensor; and the 
right panel shows the spin-separation correlation.The left 
panel of Figure 8 shows that there is a significant corre¬ 
lation of the minor axis of high-mass galaxies aligned per¬ 
pendicular to the separation vector towards other high mass 
galaxies. In the case of the reduced inertia tensor, the signif¬ 
icance of this signal is decreased, consistently with the fact 
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Figure 7. Correlations between the minor axis (rye, left for simple inertia tensor and middle panel for reduced inertia tensor) or spin 
{r)s, right panel) and separation vector as a function of comoving separation. The sample of galaxies is divided into three bins by their 
dynamical classification: Vja < 0.55, 0.55 < F/ir < 0.79 and F/cr > 0.79. Galaxies with low V/cr show a significant shape-separation 
correlation. In the left panel, the dotted lines show the alignment signal for those galaxies split by the mean mass of that population: 
logM* < 9.5 (black) and logM* > 9.5 (gray). Both low and high mass populations have a contribution, and higher mass galaxies have 
stronger alignments. 





Figure 8. Correlation between the direction of the the minor axis (rye) obtained from the simple inertia tensor (left panel) and from 
the reduced inertia tensor (middle panel) or spin of a galaxy (rya, right panel) with the separation vector to a galaxy within the same 
mass-selected sample as a function of comoving separation. High mass galaxies show a clear radial alignment trend towards other high 
mass galaxies. 


that the reduced inertia tensor yields rounder shapes for 
these galaxies. For low mass galaxies, the results are similar 
but with lower significance. For intermediate mass galaxies, 
we find a small negative correlation between the minor axis 
and the separation vector for the simple inertia tensor, and 
a small positive correlation when the reduced inertia ten¬ 
sor is used. This suggests that the reduced inertia tensor 
minor axis is correlated with the direction of the spin for 
this disc-like population, while this is not the case for the 
simple inertia tensor. We interpret this as a decreasing ten¬ 
dency of stellar particles to settle on a disc as a function of 
distance to the center of mass of the galaxy, possibly trac¬ 
ing merging structures. Indeed, satellite mergers within the 
Horizon-AGN simulation tend to redistribute their angular 
momentum significantly within the host (Welker et ah, in 
prep.). 

The overall physical picture that we get from these re¬ 
sults is the following. 

• Spin and shape alignments depend on galaxy dynamics, 
and these trends also translate into a mass dependence. 

• Spheroidal galaxies show a significant trend of radial 
alignments with respect to each other that is preserved to 
large separations and is more prominent for high mass galax¬ 
ies. The signal decreases in amplitude when the reduced iner¬ 


tia shapes are adopted, as the galaxy shapes become rounder 
and less sensitive to tidal debris in the outskirts. 

• Spin alignment trends are tangential around other 
galaxies and marginal, and seem to be better correlated with 
the reduced inertia tensor than with the simple inertia tensor 
shapes for intermediate and high Vja galaxies (correspond¬ 
ingly, also intermediate mass galaxies). 

The results presented so far only consider correlations 
between galaxies with similar properties. We also explore 
whether cross-correlations between subsets exist. Figure 9 
shows the cross-correlation of positions of spheroidal trac¬ 
ers {yja < 0.55) with the direction of the minor axis (red 
circles) and the spins (red triangles) of disc-like tracers 
{yjo > 0.55). Discs show a tendency to align their spins 
parallel to the separation vector to spheroidal galaxies (red 
dashed curve), as well as their minor axis (red solid curve). 
Notice that the direction of the spin and the minor axis of 
a disc galaxy (from the reduced inertia tensor) are very well 
correlated, and hence the red dashed and solid curves al¬ 
most lie almost on top of each other. Figure 9 also shows 
the cross-correlation of the positions of discs with the orien¬ 
tations of the minor axis (black circles) and the spin (black 
triangles) of spheroidal galaxies. These galaxies show a pref¬ 
erential elongation of their shapes towards the positions 
of discs, while they do not show significant alignment of 
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Figure 9. Cross-correlations of galaxy positions and shapes {rje) 
or spins {rjs) for the galaxy population divided in two Via bins. 
5 indicates tracers of the galaxy density field and s (m) indicates 
the sample for which the orientation is measured from the spin 
(minor axis). The black solid line shows the correlation between 
the positions of disc-like tracers of the density field (V/tr > 0.55) 
and the direction of the minor axis of the reduced inertia tensor 
of spheroidal galaxies (F/cr < 0.55). The red solid line shows the 
correlation between the positions of spheroidal galaxies and the 
direction of the minor axis of disc-like galaxies. The red dotted 
line shows the correlation between the position of discs and the 
spin of spheroidals, and the red line shows the correlation between 
the positions of spheroidals and the direction of the spin axis of 
discs. (Notice that the red lines almost overlap, showing that the 
spin and the minor axis of a disc point along the same direction.) 
These results suggest that the discs are preferentially clustered in 
the direction of the elongation of spheroidals, while they also tend 
to have their spins aligned in the direction of nearby spheroidals. 

their spins. When the simple inertia tensor shapes are used, 
the spins and shapes of discs become decorrelated and the 
significance of the shape alignment trend of discs around 
spheroidals is lost. From these results, we conclude that the 
discs are preferentially clustered in the direction of the elon¬ 
gation of spheroidals, while they also tend to have their spins 
aligned in the direction of nearby spheroidals (also traced 
by their reduced inertia tensor shapes). This picture is in 
agreement with disc galaxies living predominantly in a fila¬ 
mentary structure that follows the elongation of spheroidal 
galaxies at its knots. Furthermore, disc galaxies tend to have 
their spins aligned parallel to the direction of these fila¬ 
ments, and perpendicular to the elongation of the central 
spheroidal. Figure 10 shows a cartoon picture of alignments 
where the effect is exaggerated for visual purposes. 

We perform a series of tests to determine the signifi¬ 
cance of the measured signal of disc alignments with respect 
to spheroidal tracers. In particular, in this case we adopt a 
jackknife procedure (similar to that adopted for projected 
correlations) to estimate the uncertainty in the disc align¬ 
ment signal. The significance of the measurement is obtained 
from a test using only the diagonals of the jackknife co- 
variance matrix and without modelling the noise in this ma¬ 
trix. We find that the null hypothesis can be rejected at 
> 99.99% confidence level (C.L.) level for the alignment of 
disc spins around spheroidals (red dashed in Figure 9), and 
similarly for the alignment of the minor axes of discs around 
spheroidals (red solid fine in Figure 9). This level of signif- 


Figure 10. A cartoon picture of alignments, as interpreted from 
the results of Section 5.1. Discs live in filaments connecting el¬ 
lipticals and they tend to align their spin(s)/minor axes in the 
direction of the filament. Ellipticals tend to have their shapes (m 
represents the minor axis) aligned towards each other and to¬ 
wards the direction of the filaments. The effect of alignments is 
exaggerated for visual purposes by showing all galaxies perfectly 
aligned following the measured trends in the simulation. 

icance decreases to 92% when the simple inertia tensor is 
adopted to measure the shapes. On the contrary, the direc¬ 
tion of the spin of spheroidals is not correlated with the po¬ 
sition of discs (72% C.L. for null hypothesis rejection of the 
blue curves), but the minor axis direction of a spheroidal 
is anti-correlated with the position of discs with high sig¬ 
nificance (> 99.99% C.L. for both the reduced ans simple 
inertia tensor). Finally, we consider whether alignment sig¬ 
nals are still present when the orientation is defined by the 
direction of the major axis. We find that the disc alignment 
measurement in this case is more sensitive to the choice of 
reduced/simple inertia tensor. This result confirms that the 
simple inertia tensor is a worse tracer of the spin compared 
to the reduced inertia tensor, as the alignment signal loses 
significance in that case. On the other hand, the use of the 
simple or reduced inertia tensor does not change the sig¬ 
nificance of the alignment of spheroidals in the direction of 
discs. 

We conclude that: 

• discs show a significant tendency for tangential align¬ 
ment around over-densities traced by spheroidal galaxies, 

• spheroidals are preferentially elongated towards discs 
and other spheroidals, 

• and that spin is a good tracer of reduced inertia shapes 
for discs, but not for spheroidals. 

In the next section, we mimic observations by exploring 
projected ellipticity alignments. 

5.2 Projected correlations 

The intrinsic alignment signal is typically measured in the 
literature using the projected correlation function of galaxy 
positions and shapes (Equation 8). This quantity is read¬ 
ily accessible using shear measurements from survey galaxy 
catalogs. In this work, we also obtain the projected correla¬ 
tion functions of the density field and projected shapes, 
and wsx, where the density field is obtained from a random 
subsampling of 0.007% of the DM particles in the box. This 
subsampling guarantees a 10% convergence level in the DM 
power spectrum, which is similar to the expected level of 
convergence in determining galaxy shape (as we discussed 
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from Figure 2). Tenneti et al. (2015) used a similar approach 
with a comparable subsampling fraction. The measurement 
of UI 5 + and ws^ is advantageous in that it allows us to avoid 
modelling galaxy bias, or to make any assumption about its 
scale dependence. Also, we do not need to model peculiar 
velocities, as the DM and galaxy positions in the box are 
perfectly known. 

We measure the correlation function following 

Equation (8) for all galaxies with > 300 stellar particles in 
the simulation box and replacing the density tracers by the 
subsample of DM particles. Grid locking (see Appendix A) 
is not expected to contaminate the measurement by spher¬ 
ical symmetry. As a consequence, the S+R term is not ex¬ 
pected to contribute to this correlation and we neglect it in 
this section. Section A2 provides confirmation of these as¬ 
sumptions. We show the projected correlation functions of 
the density field and the -|- component of the shape from 
the simple inertia tensor and the reduced inertia tensor in 
the left panel of Figure 11. As expected from our results 
in Section 5.1, we find an anti-correlation between the -I- 
component of the shape and the density field that is sig¬ 
nificant at > 99.99% C.L. level for the simple inertia ten¬ 
sor. The negative sign indicates that the projected shapes of 
galaxies are elongated pointing towards other galaxies, i.e., 
that alignments are radial. We find a decreased tendency for 
alignments (47% C.L.) when using the reduced inertia ten¬ 
sor, consistently with rounder shapes and with the results 
presented in Section 5.1. The right panel of Figure 11 shows 
that the 5x correlation is consistent with null (at the ~ 65% 
C.L.). 

In Figure 12, we split the sample of galaxies with shapes 
into 5 bins of mass (left panel), V/a (middle panel) and 
u — r colour (right panel). All bins have approximately the 
same number of galaxies and the legend in each panel indi¬ 
cates the mean of the property considered for the galaxies in 
each bin. We find that galaxies in the lowest (highest) V/a 
{u — r colour) bin carry the strongest alignment signal. We 
find very similar results when splitting the galaxies by their 
g — r colour. In comparison, the split by mass results in a 
less clear identification of which galaxies are responsible for 
projected shape alignments. Figure 13 shows that there is a 
significant correlation between colour and V/cr for galaxies 
with redder colours. On the contrary, V/a and mass are not 
monotonically correlated, as shown earlier in Figure 1. 

Interestingly, the high V/a galaxies do not show any sig¬ 
nificant alignment in Figure 12. This is puzzling given the 
results presented in Figure 9. To elucidate this discrepancy, 
we compute the statistic using the same selection cuts 
as for Figure 9 and show the results in Figure 14. We find 
that, while the alignment of spheroids in the direction of the 
clustering of discs is still significant in projection, the disc 
alignment signal is diluted and consistent with null at the 
27% C.L. using the simple inertia tensor, but less so (89% 
C.L.) using the reduced inertia tensor. In the latter case, the 
signal was more significant from the orientation-separation 
correlation of Section 5.1. In projection, a tangential align¬ 
ment is only marginally present. (We remind the reader that 
gravitational lensing has the same positive sign for correla¬ 
tions measured between pairs of galaxies with large separa¬ 
tions along the line of sight). Notice that, as discussed in 
Section 5.1, a lower level of correlation is expected for the 
simple inertia tensor given that this is not a good tracer of 



Figure 13. Galaxies in the V/a vs u — r plane. There is a strong 
correlation between Vja and u — r colour at the red end of the 
colour distribution. Galaxies in Horizon-AGN have a bimodal 
distribution of colours and V/cr. The red horizontal line represent 
our fiducial cuts in V/a used in Section 5.1. 


spin alignment. Moreover, there are several reasons why the 
disc orientation correlation observed in Figure 9 can be di¬ 
luted in projection. One factor is the weighting by galaxy 
ellipticity in Equation (7). Face-on discs would carry no sig¬ 
nal in this statistic. The second reason for dilution is the fact 
that Equation (8) weights the signal in each 11 bin equally, 
while the alignment signal is expected to lose correlation 
as n increases. In comparison. Figure 9 showed the level 
of alignment as a function of three dimensional separation: 
pairs with large 11 have large r in that figure and lower cor¬ 
relation. 

Finally, we also study the cross correlation of galaxy po¬ 
sitions and galaxy shapes, considering all galaxies in the sim¬ 
ulation box. This correlation will include the effect of galaxy 
bias, compared to the DM-shape correlation. We show a 
comparison of cross correlation of -I- shapes and DM, and 
of + shapes and galaxy positions in the left panel of Fig¬ 
ure 11. We find that the galaxy-shape correlation traces the 
DM-shape correlation well within the error bars. There is a 
discrepancy in amplitude of the correlation in the first bin 
in the case of the simple inertia tensor. The excess power 
in the galaxy-shape correlation could arise from tidal debris 
on small scales or from the increased clustering compared to 
the DM, but we cannot draw firm conclusions from this com¬ 
parison. In general, the similarity between the galaxy-shape 
correlation and the DM density-shape correlation suggests 
that the bias parameter bg is not very different from unity 
for the sample of galaxies considered, and that it does not 
have any significant scale dependence. 


5.3 Modelling of alignment signal 

Early-type galaxy alignments are thought to arise due to 
the action of the tidal field. In this model, the tidal field 
contributes a small component to the projected ellipticity 
of a galaxy, given by (Catelan, Kamionkowski & Blandford 
2001) 
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Figure 11. wsj^ projected correlation function for all galaxies with > 300 stellar particles (left panel). This cut corresponds to a cut in 
mass of log(M*/M0) > 9. The same panel also presents a comparison between the correlation and the correlation between galaxy 
positions and + component of the shapes Wg+. The right panel shows the projected Sx correlation as a test for systematics. In both 
panels, results obtained with the simple inertia tensor are indicated with the blue line; while the red line corresponds to shapes obtained 
from the reduced inertia tensor. The measured points for Wg^ are arbitrarily displaced to larger radii by 5% for visual clarity in the left 
panel, and similarly for the reduced inertia tensor in the right panel. 





Figure 12. 1175 + projected correlation functions for 5 bins of mass (left panel) and V/cr (middle panel) and u — r colour (right panel). 
The legend indicates the mean of the considered property in each bin, and the 5 bins are approximately equally populated. For simplicity, 
we only show the correlations that correspond to shape measurements using the simple inertia tensor. The impact of using the reduced 
inertia tensor is shown in Figure 11. 


7(+,x) = - dl,d:,dy)S[(t)p\, (9) 

where Ci is a proportionality constant that parametrizes 
the response of the shape of a galaxy to the tidal field, cj>p 
is the Newtonian gravitational potential at the redshift of 
formation of a galaxy and 5 is a smoothing filter that acts 
to smooth the potential over the typical scale of the galactic 
halo (~ 1 Mpc)®. As a consequence, there is a correlation 
between galaxy positions and their intrinsic shapes, given 
by 

= ( 10 ) 


® Notice that we adopt a different sign convention than 
Singh, Mandelbaum & More (2015) and Tenneti et al. (2015), 
whereby alignments are negative if they are radial, and positive 
if tangential. 
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where bg is the galaxy bias, pcrit is the critical density of 
the Universe today, D{z) is the growth function (normalized 
to unity at 2 = 0) and Ps is the matter power spectrum. 
While these expressions are not strictly valid in the non¬ 
linear regime, it is customary to approximate the nonlinear 
scale alignments by replacing the linear power spectrum in 
Equation (10) (Bridle & King 2007) by its nonlinear ana¬ 
logue. Recent observational works are beginning to test this 
assumption (Singh, Mandelbaum & More 2015), but most 
constraints on the amplitude of alignments are still typically 
given in terms of the nonlinear alignment (‘NLA’) approx¬ 
imation of Bridle & King (2007). The gy. power spectrum 
is not presented because it is expected to average to null in 
projection. 


Equation (10) can be transformed to redshift space to 
give a prediction for the on-the-sky •Wg+ projected correla¬ 
tion. 
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Figure 14. Cross-correlation of galaxy positions and -h compo¬ 
nent of the shape applying the same selection as in Figure 9. The 
alignment of low V/cr galaxies in the direction of high V/a galax¬ 
ies (red) is significant after projection of the shapes and through 
the simulation box. On the contrary, the tangential alignment of 
high V/(T galaxies around low V/cr tracers of the density field 
(black) is diluted in projection, albeit the dashed line still rejects 
the null hypothesis at 89% C.L. The error bars in the case of the 
black dashed line have been artificially displaced to 5% larger 
for visual clarity. 


'Wg+{rp) = - 


Pcrit ^7) 

'r^D{z) 


POO PC 

/ dk, / 

Jq Jo 


dk± 


ki 


{ki + ki)h 


-P 5 (k, z) Sin(fc 27 rniax)i 2 (fc±rp), ( 11 ) 


where bg is the galaxy bias, is the projected radius, Ilmax 
is the line of sight distance over which the projection is 
carried out and J 2 is the second order Bessel function of 
the first kind. Notice that the Kaiser factor (Kaiser 1987; 
Singh, Mandelbaum & More 2015) is unnecessary to model 
the galaxies in the simulation, since we have access to the 
true positions of the galaxies along the line of sight. Simi¬ 
larly, we can safely neglect fingers-of-God effects arising from 
peculiar velocities in the nonlinear regime. 

Singh, Mandelbaum & More (2015) observed an excess 
of power on small scale alignments compared to the best 
fit NLA model to the LOWZ galaxy sample of the SDSS. 
In this regime, we also consider a halo model developed by 
Schneider & Bridle (2010). According to this model, the one 
halo power spectrum of galaxy positions and intrinsic shapes 
is given by 


(k, z) = -ah 


{kIPi? 

1 + (fc/p2)P3 ’ 


( 12 ) 


where an is the halo model alignment amplitude on small 
scales; pi, P 2 , and pa are fixed parameters based on fits 
by Schneider & Bridle (2010); and the projected correlation 
function is then given by 


wlX =-bg j '^^k±Pl'^^i{k^,z)Jo{ki_rj,). (13) 


While the alignments of disc galaxy shapes have not 
been detected in observations, and despite the fact that 
the tidal alignment model is not expected to describe 


their alignments (Catelan, Kamionkowski & Blandford 
2001; Hirata & Seljak 2004), the available constraints are 
typically phrased in terms of the NLA model as well (e.g., 
Mandelbaum et al. 2011). 

We now place quantitative constraints on the goodness- 
of-fit of the intrinsic alignment models. We build a model 
template of the intrinsic alignment signal, Wg_f_(rp), at the 
positions for which we measure this correlation in the sim¬ 
ulation. We model the nonlinear matter power spectrum 
using the CAMB software (Hewlett et al. 2012) with the 
HALOFIT correction (Smith et al. 2003). 

Using the diagonals of the jackknife covariance, we de¬ 
termine the from summing over all radial bins as follows 

2 _ (w6-r(rp) — 

^ ^ Var[w 5 +(rp)] 

where p represents the alignment model template to be fit 
and we look for the minimum x^ by varying the parameters 
of the fit; and analogously for rcg+, in which case only joint 
constraints of the product of bg and the alignment ampli¬ 
tude are obtained. We emphasize that we are not modelling 
covariances or noise in this matrix when performing the fits. 

Notice that while the results presented in the previ¬ 
ous section do not guarantee that disc-like galaxies do not 
contribute to the alignment signal, the correlation functions 
shown in Section 5.1 suggest that they counter-act the radial 
alignment of ellipticals due to their tendency to orient tan¬ 
gentially around over-densities. For this reason, and given 
that we do not have sufficient constraining power to bin 
the galaxy sample into different populations, we fit different 
models to Wg+ and shown in Figure 11, which include 
the contribution of all galaxies with > IO^Mq stellar masses 
in the simulation. 

We focus on fits to the signal obtained using the simple 
inertia tensor, as in the case of the reduced inertia tensor the 
results are consistent with null. A power-law fit as a func¬ 
tion of projected radius (rCg_,_ = bgAiVp) yields the following 
constraints on power-law amplitude and power-law index: 
bgAi = —0.13 ± 0.03 and /3 = —0.75 ± 0.25, respectively. 
The error bars quoted correspond to 68 % C.L when holding 
one of the parameters fixed at the best fit value correspond¬ 
ing to the minimum x^ {bgAi = —0.17 and /3 = —1.0). We 
find that the linear alignment model is a very poor fit due to 
the lack of power on small scales compared to the measured 
signal. Similarly, the NLA model also underestimates the 
correlation of positions and shapes on small scales^. In this 
case, the best fit model is shown in Figure 15. The best fit 
NLA model for Wg+ in the case of the simple inertia tensor is 
represented by the black solid line. In gray, we show the best 
fit NLA model for ' 015 +, also for the simple inertia tensor, 
which is comparable to the black line. The black dotted line 
shows the best fit to ws+ in the case of the reduced inertia 
tensor, for which the best fit amplitude is consistent with 
null at 68 % C.L. 

Given that the best fit NLA model underestimates the 
alignment signal on small scales, we consider fitting a sum 


In this particular work, we do not apply any smoothing to 
the tidal field. A smoothing filter would suppress power on small 
scales, worsening the comparison of the NLA model to the mea¬ 
sured signal. 
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of the NLA model template (across all scales measured) and 
the halo model (only at Vp < Q.8/h) with the parameters of 
Schneider & Bridle (2010). Notice that the physical inter¬ 
pretation of these results is not straightforward and should 
only be considered as phenomenological. The constraints on 
the parameters for the NLA and halo model combination 
are: bgan — 0.271q;ii, hgAi = O.SSlg gg. The in this case 
is comparable to that of the power-law fit. If we use more 
and smaller (I//4 a side) jackknife regions, the error bars 
increase slightly yielding bgAi = 0 . 60 lg; 4 j. 

Given the limited simulation volume and the corre¬ 
sponding (cosmic) large error bars, we cannot obtain mean¬ 
ingful constraints on mass or luminosity dependence of the 
signal. However, we note that, as discussed in Section 5.2, 
the mass dependence is not monotonic. 


6 DISCUSSION 

Theoretical models of intrinsic alignments have in¬ 
deed suggested that the population of galaxies 
subject to this mechanism can be split in two 
(Catelan, Kamionkowski & Blandford 2001). Disc-like sys¬ 
tems interact with the large-scale structure through torques 
to their angular momentum vector, while spheroidals, which 
do not have significant angular momentum, would tend to 
orient their major axes pointing towards over-densities. In 
this work, we have used U/cr as a proxy of galaxy morphol¬ 
ogy, and we have indeed confirmed the existence of two 
different alignment mechanisms in play for spheroidals (low 
V/a) and disc-like (high U/n) galaxies. While Vja is not 
easily accessible to upcoming gravitational lensing surveys, 
fortunately there is a strong correlation between u — r and 
U/(T (Figure 13) that can help identify the two populations 
of galaxies and their alignments. Low U/cr galaxies tend 
to have redder colours, while high U/cr galaxies tend to be 
bluer. 

The alignments of blue galaxy shapes have so far been 
consistent with null from observations (Mandelbaum et al. 
2011; Heymans et al. 2013), albeit with large error bars 
that still allow for a significant level of contamina¬ 
tion from blue galaxy alignments to cosmological observ¬ 
ables in current and future surveys (Chisari et al. 2015; 
Krause, Eifler & Blazek 2015). We similarly find that the 
projected correlation function of blue galaxy positions and 
shapes is consistent with null in Horizon-AGN. Neverthe¬ 
less, while blue galaxies do not align around each other, they 
tend to align tangentially around red galaxies in three di¬ 
mensional space, and this signal is washed out in projection 
within our error bars. We emphasize that this does not im¬ 
ply that blue galaxy alignments can be neglected for future 
surveys. Codis et al. (2015) found that spin-spin alignments 
could potentially translate into worrisome levels of contami¬ 
nation for future lensing surveys, particularly for blue galax¬ 
ies. Their results were based on a “spin-gives-ellipticity” pre¬ 
scription for translating spins into shapes for disc-like galax¬ 
ies. We found that the validity of this assumption depends on 
the method used to determine ellipticity, with better agree¬ 
ment for the reduced inertia tensor case. The statistics of 
spin/shape alignments of disc galaxies in three dimensions 
from the simulation can be a useful tool to constrain disc 


alignment models and determine the level of contamination 
to future surveys. 

Several theoretical predictions have been made for the 
spin alignments of disc-like galaxies. In simulations, discs 
tend to form with their spin correlated with the direction 
of filaments (Bailin & Steinmetz 2005; Aragon-Calvo et al. 
2007; Hahn et al. 2007; Sousbie et al. 2008; Zhang et al. 
2009; Codis et al. 2012; Libeskind et al. 2012; Dubois et al. 
2014) and the vorticity of the density field around them 
(Laigle et al. 2015). This is consistent with the scenario of 
spin acquisition by tidal torquing biased by the large scale 
structure filaments. Extended tidal torque theory predicts 
that halos with mass below 5 x 10^^ Mq at redshift 2 = 0 
tend to form with their spin pointing along the direction 
of filaments, the center of which are represented by saddle 
points of the density field (Codis, Pichon & Pogosyan 2015). 
In the plane perpendicular to the filament and containing 
the saddle point, halos are predicted to show this prefer¬ 
ential orientation. Away from this plane and closer to the 
nodes of the filament where higher mass halos reside, spins 
flip direction, becoming perpendicular to the filamentary 
axis. Hence, if galactic spins are correlated with halo spins, 
we should expect low-mass galaxies to have a spin aligned 
with the filamentary axis (as found in Dubois et al. (2014) 
for Horizon-AGN galaxies) and therefore also aligned® with 
the separation vector (since we expect a higher number 
of galaxies to inhabit within the filament). In contrast, 
Codis, Pichon & Pogosyan (2015) predict that very massive 
galaxies should have a spin perpendicular to the filament. 
However, massive galaxies tend to be supported by random 
motions of their stars, and hence the definition of the spin 
direction becomes more noisy. The signal predicted by the¬ 
ory could also be diluted due to the fact that not all low 
mass galaxies reside in filaments. We found that the spin 
and the minor axes of massive galaxies have less tendency 
to be aligned than for their lower mass counterparts. 

Works by Tenneti et al. (2014, 2015); Velliscig et al. 
(2015) have also studied the intrinsic alignments of galax¬ 
ies in hydrodynamical simulations. Tenneti et al. (2015) find 
that orientation-position correlations (their ‘ED’ correla¬ 
tion) with the density field have similar strengths and signs 
for blue and red galaxies when correlated with the density 
field. On the contrary, we find that blue galaxies tend to have 
tangential alignments around the locations of red galaxies 
(but not around each other). In projection, both works find 
the blue galaxy alignments to be suppressed. Those authors 
did not study spin alignments, which in this work have been 
shown to be significant for disc-like galaxies and to be con¬ 
nected to their shape alignments. 

Observational results have shown there is a 
strong trend for radial alignment of the ellipticities of 
red galaxies (Mandelbaum et al. 2006a; Hirata et al. 
2007; Okumura, Jing & Li 2009; Joachimi et al. 2011; 
Heymans et al. 2013; Singh, Mandelbaum & More 2015), 
which increases for higher mass galaxies. This trend is 
also seen in Horizon-AGN for the shapes of spheroidals. 
Tenneti et al. (2015) presented NLA and power-law fits to 
ws+ from their simulation as a function of redshift and 


° or anti aligned, should the galaxies belong to octants of opposite 
polarity, see Godis, Pichon &; Pogosyan (2015), section 6.2.3. 
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Figure 15. Alignment model fits for and rOg+ for all galaxies with > 300 stellar particles. Results obtained with the simple inertia 
tensor are indicated with the red lines; while the blue line corresponds to shapes obtained from the reduced inertia tensor for 105 + only. 
We show fits from the NLA model to ws+ (gray solid for simple inertia tensor; black dotted for reduced inertia tensor) and to Wg+ (black 
solid for the simple inertia tensor). The measured points for and simple inertia tensor shapes are arbitrarily displaced to larger 

radii by 5% for visual clarity. The NLA model significantly underestimates the power in alignments at small separations (< 1 Mpc//i) in 
qualitative agreement with observations by Singh, Mandelbaum &; More (2015). 

luminosity of the galaxy sample using iterative reduced 
inertia tensor for the galaxy shapes. We find that their 
fits reproduce the alignment signal of luminous galaxies 
{Mr < —22.6) presented in this work. We show this agree¬ 
ment in Figure 16. The gray curves show their power-law fits 
on small scales, and their NLA model fits across all scales 
(see their Table 1). Given that Tenneti et al. (2015) match 
the observed alignments of LRGs with their simulation 
data, it is expected that Horizon-AGN will equally match 
observations if the redshift dependence of the signal is 
similar to that found in that work. However, we find that 
their fits significantly overestimate the alignment amplitude 
of the whole galaxy sample presented in Figure 15, and 
this could indicate a steeper luminosity dependence of 
alignments in Horizon-AGN. Joachimi et al. (2011) indeed 
found a steeper luminosity dependence of observed LRG 
alignments compared to Tenneti et al. (2015). However, 
such a steep luminosity dependence is not sufficient to 
reproduce the alignment amplitude of Figure 15. It is 
possible that blue galaxy alignments are suppressing our 
results in that figure as well. 

Godis et al. (2014) found no alignments for red galax¬ 
ies using Horizon-AGN at z = 1.2. In that work, spin was 
used as a proxy for galaxy shape. Our results are consistent 
with theirs. We have shown that for red galaxies, the spin 
alignment signal is very different from the shape alignment 
signal, as these galaxies do not carry significant angular mo¬ 
mentum. 

Velliscig et al. (2015) studied the galaxy-halo misalign¬ 
ment comparing the shapes of stars and hot gas to that of the 



Figure 16. Measurement of ws+ from Horizon-AGN for the 
most luminous galaxies, with absolute r-band magnitudes Mr < 
—22.6 (~ 800 galaxies). The blue curve corresponds to shapes 
measured using the reduced inertia tensor and the red curve, us¬ 
ing the simple inertia tensor. We also show alignment model fits 
from Tenneti et al. (2015) (gray), which are in good agreement 
with our results. Tenneti et al. (2015) fit a power-law on small 
scales (0.1 — 1 Mpc/h, shown here in the range 0.1 — 2 Mpc/h) and 
the NLA model on large scales (6 — 25 Mpc//i, extrapolated here 
over all scales shown). While the agreement is good at high lumi¬ 
nosity, we find that the fits by Tenneti et al. (2015) significantly 
overpredict the alignment signal for the whole sample shown in 
Figure 15. 


O 
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underlying DM halo using the EAGLE smoothed-particle- 
hydrodynamics simulation. They found that the alignment 
of the stellar component with the entire DM halo increases 
as a function of distance from the center to the subhalo and 
as a function of halo mass. This is in qualitative agreement 
with our finding that the shape alignment signal is reduced 
when using the reduced inertia tensor, which puts more 
weights towards the inner regions of galaxy. However, they 
also find that misalignment angles between the stellar and 
the DM components are larger for early-type than for late- 
type galaxies using the simple inertia tensor. This result will 
require further comparison, as the DM is expected to have 
stronger alignments than the baryons (Okumura, Jing & Li 
2009) and given that we find a stronger shape alignment 
signal in Horizon-AGN for early-types than for late-types. 


7 CONCLUSIONS 

We have studied the alignments of galaxies, as traced by 
their stellar particles, using the Horizon-AGN simulation. 

The main result of this paper is the clear identifica¬ 
tion of two different alignment mechanisms for disc-like 
galaxies and spheroidals, in qualitative agreement with the¬ 
oretical expectations (Catelan, Kamionkowski & Blandford 
2001; Codis, Pichon & Pogosyan 2015). This is the first 
time that these two mechanisms are clearly separated in a 
hydrodynamical cosmological simulation. In contrast, pre¬ 
vious work by Tenneti et al. (2014) was unable to dis¬ 
tinguish between red/blue galaxy orientation-separation 
alignments. This is likely a consequence of the different 
methods used to solve for the hydrodynamics (AMR in 
Horizon-AGN compared to smoothed particle hydrodynam¬ 
ics in MassiveBlack II), which result in different galaxy prop¬ 
erties and their evolution with redshift. 

We also reached the following conclusions: 

• High mass (low V/a) galaxies are elongated pointing 
towards other galaxies. This trend is preserved when pro¬ 
jected correlations of the density field and galaxy shapes 
are considered. 

• There is a preferential tangential orientation of disc¬ 
like galaxies around spheroidals. This trend is diluted in 
projection, possibly due to the equal weighting of the dif¬ 
ferent H bins in Equation (8) and the reduced contribution 
of face-on discs to Equation (7). This suggests that in order 
to extract the maximum possible information from intrinsic 
alignments in simulations, it could be beneficial to avoid per¬ 
forming projections along the line of sight, and rather access 
the full three dimensional information on alignments pro¬ 
vided by the simulations. This would increase the signal-to- 
noise ratio in intrinsic alignments constraints from simula¬ 
tions and allow for more accurate forecasts of intrinsic align¬ 
ment contamination to weak leasing in future surveys, such 
as Euclid^ (Laureijs et al. 2011), the Large Synoptic Survey 
Telescope^(Ivezic et al. 2008) and WFIRST® (Green et al. 
2011 ). 


® http://sci.esa.iiit/euclid 
^ http://www.lsst.org 
® http://wfirst.gsfc.nasa.gov/ 


• We are able to describe Wg+ across all scales probed 
and for all galaxies in the simulation with > 10® Mq using 
a power-law in the case where shapes are obtained with the 
simple inertia tensor. The NLA and LA model tend to under¬ 
estimate the power on small scales; a conclusion also reached 
by Singh, Mandelbaum & More (2015) using low redshift 
observations of LRGs. Fits to in the reduced inertia 
tensor case result in an alignment amplitude that is con¬ 
sistent with null for our complete sample of galaxies. This 
does not imply that alignments are not potential contam¬ 
inants to weak leasing measurements. It will be necessary 
to match the shape measurement and galaxy selection done 
in observations to make more quantitative assessments of 
contamination from alignment to future surveys. The align¬ 
ments of the most luminous galaxies in Horizon-AGN are, 
in fact, in agreement with work by Tenneti et al. (2015) and 
the alignments of bright LRGs in SDSS. 

• Galactic kinematics, as quantified by V/a, is a good 
proxy for the level of alignment. We also find that, given the 
existing correlation between V/a and u — r colour, the latter 
can also be used as proxy to separate galaxy populations 
with different sensitivity to alignment. On the contrary, the 
amplitude of the intrinsic alignment signal is not monotonic 
with stellar mass, which we interpret as a consequence of 
the wide distribution of stellar masses at low V/a. 

• We emphasise that correlations of spins (or shapes) and 
separation are not contaminated by grid-locking because 
this effect averages to null for position-shape correlations 
(as shown in Appendix A). This is not necessarily true for 
two-point auto-correlations of spins and/or shapes, and so 
we refrain from giving these a physical interpretation in the 
main body of the manuscript, although we present them in 
Appendix B for completeness. 

• The choice of shape estimator can have a large impact 
on predictions for intrinsic alignments from cosmological 
simulations. To make accurate predictions of intrinsic align¬ 
ment impact on weak leasing surveys, we expect that it will 
be necessary to create mock images of simulated galaxies in 
a manner that takes into account photometric depth, noise 
and convolution by the point-spread function. We defer this 
and a study of redshift evolution and selection cuts on the 
galaxy sample to future work. 
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APPENDIX A: GRID-LOCKING 

Adaptive-mesh-refinement simulations are thought to be 
subject to grid-locking systematics, whereby the galaxy 
spins become aligned with the directions of the grid. We 
have verified that this effect is indeed present in the 
Horizon-AGN simulation, but we demonstrate in this ap¬ 
pendix that it does not affect position-shape or position-spin 
correlations. 


A1 Correlations with the box 

We consider possible systematics that might arise from cor¬ 
relations of the galaxy orientations with preferential direc¬ 
tions of the simulation box. Figure A1 shows the excess frac¬ 
tion of galaxies as a function of | cos(s ■ {x, y,z})|, where 
s is the galaxy spin and we average over x, y, z. There is 
a clear excess of galaxies at directions perpendicular and 
parallel/anti-parallel to the grid, and a decrement at inter¬ 
mediate angles. A similar behaviour is observed for the direc¬ 
tion minor axes in the same figure. The higher mass galax¬ 
ies are less affected, but the trend is not monotonic with 
mass. We have also considered selecting the galaxy sample 
by median stellar age, resolution and colour. While mass 
and median stellar age cuts can help remove the system¬ 
atics, it is not possible to define clean sample without re¬ 
moving the vast majority of the galaxies. Moreover, it is not 
clear how these selection cuts affect the intrinsic alignment 
measurement and its comparison to current observational 
constraints. For these reasons, we decide to avoid placing 
cuts on the galaxy sample to reduce the grid-locking effects. 
Instead, we show in the next section that the correlation be¬ 
tween galaxy shapes, spins and the box have no impact on 
the W 5 + statistic presented in Section 4. 

We have also considered the grid-locking signal of 
shapes obtained through the reduced inertia tensor. In this 
case, we find that the grid-locking signal is smaller than for 
Figure Al. 


A2 Gontamination to intrinsic alignments 

Grid-locking creates correlations of the spins and shapes 
of galaxies with {x, y, z}. However, due to the periodic 
boundary conditions of the box, we expect that the grid¬ 
locking will average to null around any arbitrary point 
within the box. We test this hypothesis by randomizing the 
position of galaxies while computing the relative orienta¬ 
tion of their spins and shapes around these random posi¬ 
tions. Figure A2 shows an example correlation for shapes 
defined from the orientation of the minor axis and galax¬ 
ies in three mass bins. We do not find significant correla¬ 
tions when positions are randomized. Notice that this test 
is sufficient only under the assumption that the large-scale 
structure is not itself locked to the grid. This assumption 
was shown to be valid for the dark matter large-scale struc¬ 
ture in Horizon-AGN (Dubois et al. 2014). We further test 
this assumption for galaxies in Horizon-AGN in Figure A3, 
where we show the excess probability of spin alignments (top 
panel) and separation vectors alignments (bottom panel) 
with the grid for pairs of galaxies separated by less than 25 
Mpc/h in the simulation. Spins are clearly correlated with 



Figure A2. Gorrelation between minor axis orientation and 
separation vector as a function of separation using random¬ 
ized positions. Black lines correspond to low mass galaxies 
(10® Mq < M* < 10® ® Mq); red lines, to intermediate mass 
galaxies (10® ® Mq < M* < 10^®'® Mq); and blue lines, to high 
mass galaxies (M* > 10^®'® Mq). 

the box axes, as demonstrated in Fig. Al, while this is not 
the case of the separation vectors. 

We apply a similar procedure to confirm that grid¬ 
locking does not contribute to ws+. We obtain the pro¬ 
jected correlation function of the random positions and the 
-b (inr+) and x components of the shapes. Both of these cor¬ 
relations are consistent with null. This confirms that spher¬ 
ical averages remove the effect of grid-locking on position- 
shape correlations. Figure A4 shows Wr+ for five mass bins 
as a function of Xp. 

Notice that this procedure is not applicable to two-point 
auto-correlations of the spin or the minor axes. In other 
words, the impact of grid-locking on auto-correlations can¬ 
not be fully quantified by randomizing galaxy positions. This 
procedure would only measure the impact of uniform corre¬ 
lations of galaxy shapes or spins across {x,y,z}, but two 
galaxies clustered together could be increasingly affected by 
grid-locking. We present auto-correlations in Appendix B, 
along with possible interpretations of the signals as a result 
of physical alignments or grid-locking. 


APPENDIX B: AUTO-GORRELATIONS 

We have so far shown that the spins and shapes of galax¬ 
ies tend to align with the grid, and that these grid-locking 
effects do not affect our measurements of position-shape 
(or position-spin) correlations. Unfortunately, the symme¬ 
try arguments that apply to position-shape correlations no 
longer hold when auto-correlations are considered. As a con¬ 
sequence, it is harder to separate the physical alignment 
signal and the grid-locking from auto-correlation of galax¬ 
ies spins or shapes. In this section, we present those auto¬ 
correlations and we give the two possible interpretations of 
the signal. 

The spin-spin (SS) correlation function measures the 
relative orientation of the angular momenta of two galaxies 
separated by a comoving distance vector r, 

r?s(r) = (|s(x) ■ s(x-f r)|^) - 1/3, (Bl) 

and the ellipticity-ellipticity (EE) correlation function is 
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Figure Al. Excess fraction of galaxies with respect to a uniform distribution as a function of the angle between the spin/minor axis 
with the directions of the box (averaged over j = {x^y^z}). In the left panels, the galaxy sample is split by mass in the following bins: 
10® Mq < M* < 10®'^ Mq (black), 10®'^ Mq < M* < 10^® ® Mq (red) and M* > 10^® ® Mq (blue). In the left panel, the galaxy sample 
is split by V/cr: Vju < 0.55 (black), 0.55 < Vju < 0.79 (red) and V/cr > 0.79 (blue). There is an excess of galaxies in directions parallel 
and perpendicular to the grid axes, which is less significant for high mass and low V/cr galaxies. This effect averages to null for the 
angular and projected correlations of spin/minor axis and separation, but it could give rise to a non-negligible two-point correlation. 



Figure A3. Distribution of spins (top panel) and separation vec¬ 
tors of galaxy pairs separated by less than 25 Mpc/h (bottom 
panel) on the sphere. The excess probability ^ defined so that the 
PDF reads P(cos0,0) = (1 -|- (^)/47r is colour-coded from dark 
blue (-1) to dark red (+1). If the spins are clearly correlated with 
the box axes, this is not the case of the separation vectors for 
which the magnitude of the excess probability is smaller (lighter 
colours) and the pattern does not seem to be correlated with the 
grid. 


given by 

T?is(r) = (|e(x) ■e(x + r)|^) - 1/3, (B2) 

where e is the direction of the minor axis. Similarly to ^g+, 
defined in Equation (6), we can define the auto-correlation 
function of -|- components of the ellipticity 

^++(^P,n) = ^, (B3) 

and its projection along the line-of-sight, rc++; and analo¬ 
gously for ^xx and Wxx- 

Figure B1 shows the relative orientation of the spins of 



Figure A4. Projected correlation function between random posi¬ 
tions and -|- component of the ellipticity. This result is consistent 
with null, which suggests that any systematics coming from grid¬ 
locking is not affecting our 5+ or (5x measurements. The different 
colours represent different bins in stellar mass and the points are 
slightly displaced to higher r-p for each mass bin for visual clarity. 
The mean stellar mass in each bin is indicated in the legend. 


two galaxies (right panel) separated by a comoving distance 
of r. The correlation is significant for low and intermediate 
Vja galaxies. For pair separations < 10 h~^ Mpc, the spins 
tend to align with each other and the alignment decreases on 
larger scales. This trend is monotonic with V/a. The results 
for the auto-correlation between minor axes or between spins 
are very similar when minor axes are computed using the re¬ 
duced inertia tensor (middle panel of Figure Bl). The signal 
slightly decreases when the simple inertia tensor is used (left 
panel). These three dimensional auto-correlations of galaxy 
orientations seem to be dominated by the contribution of 
disc-like galaxies. In Section 5.1, we saw that disc galaxies 
tend orient their spin (also reduced minor axis) tangentially 
around over-densities. Figure Bl is qualitatively compatible 
with those results. Similarly, we find no auto-correlation be- 
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tween the shapes of low V/a galaxies in the left panel of 
Figure Bl. 

We show the orientation-separation correlations as a 
function of mass in Figure B2. Intermediate and low mass 
galaxies tend to have a stronger auto-correlation than 
high mass galaxies. The signal is stronger for interme¬ 
diate mass galaxies, which is expected both from grid¬ 
locking and from the physical signal due to the contribu¬ 
tion of the high Vja population. Similar qualitative results 
were obtained by Codis et al. (2015) at z = 1.2 using the 
Horizon-AGN simulation with the same mass selection (see 
their Figure 8 for a direct comparison). The trend with mass 
is similar, but the amplitude of the signal is higher in our 
case, suggesting that the spin auto-correlation increases at 
lower redshift and at fixed mass. Codis et al. (2015) also 
found that separating the galaxy population hj u — r colour 
into red and blue galaxies resulted in a significant spin align¬ 
ment signal for blue galaxies, while no spin alignment was 
found for red galaxies. This is also in agreement with our re¬ 
sults, given the strong correlation between colour and V/a. 

The projected auto-correlation functions of galaxy 
shapes in three mass bins and three Vja are shown in Fig¬ 
ure B3 decomposed into ill-modes (-1—b) and B-modes (x x). 
Qualitatively, we observe that the projected auto-correlation 
of shapes is more significant for intermediate mass galax¬ 
ies and high V/a. This is consequence of a selection effect: 
galaxies with intermediate masses also tend to have higher 
V/a in Figure 1. We also find that the amplitude of the 
auto-correlation decreases when the reduced inertia tensor 
is used. We attribute this to the fact that w++ and Wxx 
carry a double weighting by galaxy ellipticities, which are 
also more round in the reduced case. 

Both B-modes and B-modes are present with similar 
amplitudes and scale-dependence in the top middle and right 
bottom panels of Figure B3. Using the diagonals of the jack¬ 
knife covariance, we find that the signal is different from null 
at the 89% C.L. (85% C.L.) for the B-modes (B-modes) of 
intermediate mass galaxies when using the simple inertia 
tensor; and higher (99% C.L.) using the reduced inertia ten¬ 
sor. For high V/a galaxies, the signal is different from null 
at the 79% C.L. for B-modes and 81% for B-modes using 
the simple inertia tensor; and higher (> 93% C.L.) using the 
reduced inertia tensor. All other bins of mass and V/a have 
auto-correlations generally consistent with being null at the 
2a level. 

Overall, we find that auto-correlations of galaxy shapes 
or spins are dominated by the contribution of disc-like galax¬ 
ies. We find no contribution from spheroidals to these auto¬ 
correlations. Moreover, we find equal contribution of B- 
modes and B-modes in the projected correlation functions. 
If the signal is to be interpreted as a physical signal, the pres¬ 
ence of an auto-correlation of disc-like tracers is expected, 
and the trend with mass/dynamics agrees with that found 
in Section 5.1. 

On the other hand, the absence of an auto-correlation 
for the shapes of spheroidals seems inconsistent with the 
presence of auto-correlations for discs. This could suggest 
that at least part of this signal is produced by grid-locking. 
B-mode and B-mode auto-correlations could also be pro¬ 
duced by both a physical signal or by grid-locking. In the 
latter case, B-modes would be a consequence of preferen¬ 
tial alignments of the simulated galaxies with the diago¬ 


nals of the grid. Ideally, two-point correlations such as those 
presented in Figures Bl and B3 could potentially be used 
to constrain a simulation-dependent model of grid-locking. 
However, there is not sufficient evidence that current intrin¬ 
sic alignment models predict the correct relation between 
density-shape and shape-shape correlations, especially on 
such small scales as probed in this work. Hence separat¬ 
ing the grid-locking from the physical alignment signal in 
auto-correlations remains a difficult task for adaptive-mesh- 
refinement codes like ramses. 
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Figure Bl. Correlations between minor axes obtained from the simple inertia tensor (left panel), from the reduced inertia tensor (middle 
panel) and between spins (right panel) as a function of comoving separation. Black lines correspond to low V/a galaxies {V/a < 0.55); 
red lines, to intermediate V/a galaxies (0.55 < V/a < 0.79); and blue lines, to high V/a galaxies {V/a > 0.79). 





Figure B2. Correlations between minor axes obtained from the simple inertia tensor (left panel), from the reduced inertia tensor 
(middle panel) and between spins (right panel) as a function of comoving separation. Black lines correspond to low mass galaxies 
(10® Mq < M* < 10®'^ M©); red lines, to intermediate mass galaxies (10®'^ M© < M* < 10^® ® Mq); and blue lines, to high mass 
galaxies (M+ > 10^®'® M 0 ). 


Low mass 



Intermediate mass 







Figure B3. Projected auto-correlation functions of -h and x components of the shape as a function of projected radius in three mass 
bins (upper panels) and three V/a bins (lower panels). In all panels, the H—h correlation is shown in black for shapes measured using 
the simple inertia tensor; correspondingly, the xx correlation is shown in red. All panels show the impact of using the reduced inertia 
tensor (blue lines for and green lines for xx). All xx correlations are arbitrarily displaced to 5% higher for visual clarity. 
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